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Abstract 

Our Universe is ruled by quantum mechanics and its extension Quantum Field Theory (QFT). 
However, the explanations for a number of cosmological phenomena such as inflation, dark energy, 
symmetry breakings, and phase transitions need the presence of classical scalar fields. Although the 
process of condensation of scalar fields in the lab is fairly well understood, the extension of results to 
a cosmological context is not trivial. Here we investigate the formation of a condensate - a classical 
scalar field - after reheating of the Universe. We assume a light quantum scalar field produced 
' by the decay of a heavy particle, which for simplicity is assumed to be another scalar. We show 

CSJ , that during radiation domination epoch under certain conditions, the decay of the heavy particle 

\^ ' alone is sufficient for the production of a condensate. This process is very similar to preheating 

, ^ , - the exponential particle production at the end of inflation. During matter domination epoch 

when the expansion of the Universe is faster, the decay alone can not keep the growing trend of 
the field and the amplitude of the condensate decreases rapidly, unless there is a self interaction. 
This issue is particularly important for dark energy. We show that quantum corrections of the 
self-interaction play a crucial role in this process. Notably, they induce an effective action which 
includes inverse power-law terms, and therefore can lead to a tracking behaviour even when the 
classical self-interaction is a simple power-law of order 3 or 4. This removes the necessity of having 
nonrenormalisable terms in the Lagrangian. If dark energy is the condensate of a quantum scalar 
field, these results show that its presence is deeply related to the action of quantum physics at 
(— I ■ largest observable scales. 



in 



> ■ 1 Introduction 
0^ 

0> ■ Observations of phenomena such as superconductivity and super fluidity in condense matter indi- 

! cates that quantum particles can collectively behave like a classical self-interacting scalar field. The 

' potential energy of this interaction plays an important role in breaking global and/or local (gauge) 

■ symmetries which usually are followed by a phase transition. The same phenomena is assumed to hap- 



o 



pen at fundamental level in particle physics where usually a quantum scalar field, e.g. Higgs boson is 
responsible for dynamical mass generation. Other phenomena, mostly cosmological such as inflation, 
leptogenesis, and many of candidate models for dark energy are based on the existence of a classical 
scalar field which is usually related to a fundamental quantum scalar field because the physics of the 



X 

H , 

5r , Universe and its content in its most elementary level is quantic. 

A classical field is more than just classical behaviour of a large number of scalar particles. In a 
quantum system particles can be in superposition states i.e. quantum mechanically correlated to each 
others. Decoherence which is generated by interaction of each particle or field with its environment 
remove the quantum superposition and correlation between quantum states, but this does not mean 
that after decoherence of scalar particles, they behave collectively like a classical field. The following 
simple example can demonstrate this fact: 

Consider a closed system consisting of a macroscopic amount of unstable massive scalar particles 
which decay to a pair of light scalar particles with a global SU (2) symmetry and a very weak coupling 
with each other. If the unstable particle is a singlet of this symmetry, the remnant particles are 
entangled by their SU{2) state. After a time much larger than the lifetime of the massive particle, 
the system consists of a relativistic gas of pair entangled particles. If a detector measures this SU{2) 
charge without significant modification of their kinetic energy, the entanglement of pairs will break 
i.e. the system decoheres and becomes a relativistic gas. The equation of state of a relativistic ideal 
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gas is Wrei = P/p w 1/3 where the pressure P and density p are defined as the expectation value of 
some operators acting on the Fock space of the system. By contrast, a classical scalar field (p{x) is a 
C-number and its density p^, pressure P^^, and kinetic energy are defined as: 



p^ = K^ + V{ip) 
P^ = K^- Viip) 



(1) 
(2) 



(3) 



where V{(p) is a potential presenting the self-interaction of the field <f{x). When it is much smaller 
than kinetic energy , we obtain » p^, and if V{ip) » K^,, — Pip- Therefore in general, a 

relativistic gas and a scalar field do not share the same equation of state, and the proof of decoherence 
in a system is not enough when a classical scalar field is needed to explain a physical phenomenon. 

Historically, the concept of a classical scalar field was first appeared in the context of scalar-tensor - 
Brans-Dicke - gravity theories (see [T] for a historical review) . In these models the scalar field presents 
dilaton, the generator of conformal symmetry. Therefore it had a purely geometric nature. It was 
only later when people tried to quantize Einstein and other gravity models that this field got a particle 
interpretation. The discovery of Higgs mechanism and other phenomena in condense matter in which 
scalar fields are present, encouraged this interpretation. More recently scalar field are found to be a 
principle ingredient in supersymmetric and superstring theories. In the classical limit quantum scalar 
fields are usually identified with classical fields and their differences are overlooked. 

When a classical system is quantized, according to canonical quantization procedure, classical observ- 
ables are replaced by operators acting on a Hilbert or Fock space, respectively for a single particle and 
for a multi-particle quantum system. The expectation values of these operators are the outcome of 
measurements. Therefore, it is natural to define the classical observable related to a quantum scalar 
field as its expectation value: 



where |^') is the state of the quantum system i.e. an element of the Fock space of the system. In 
analogy with particles in the ground state in quantum mechanics, the classical field ip{x) is also called a 
condensate. In fact a coherent state consisting of superposition of particles in the ground state behaves 
like a classical field i.e. (^'|$(a;)|^) ^ [2]. This is an ideal and exceptional case in which the number 
of particles in the system is infinite. Nonetheless, in the cosmological context where the number of 
particles is very large it can be a good approximation. Thus, later in this work we use this state to 
calculate the evolution of a condensate in an expanding universe. Using canonical representation, it is 
easy to see that for systems with a limited number of scalar particles = 0. But in presence of 

an interaction, even after renormalization, a finite term can survive [13] to play the role of a classical 
field (condensate) according to the definition In fact $ can be considered to be dressed and its 
expectation value even on the vacuum can be non-zero. Equivalently, ^ can be considered as a free 
field. In this case |\I') must include infinite number of interacting particles. In both interpretations 
the presence of an interaction is a necessary condition for the condensation of a finite system [3] (see 
also Ref. [H [23] for a review) . 

Although classical scalar fields play crucial roles in the modeling of many phenomena particle physics 
and cosmology, their existence is usually considered as granted and the efforts are concentrated on 
the relevant potentials, solutions of their dynamic equations, and quantization of small fluctuations 
around the classical background fields. For instance in the context of inflation and reheating of the 
Universe, fluctuations of inflaton are quantized around the uniform and classical background which 
is responsible for the exponential expansion of the Universe (see e.g. Ref. [5] for a review). Both in 
inflation and in ultra-cold matter the presence of a condensate is apriori justiflable. If the entropy of 
the system e.g. Universe before inflation was very small and inflatons were the dominant content, most 
of them had to be in their ground state - the zero mode - and therefore according to Ref. |2] behaved 
as a condensate (see also Appendix [Bj for a more general description of a condensate state) . In other 
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contexts such as in the reheating era, and in cosmological and lab phase transitions the entropy is not 
always small. Therefore, as the above example showed, in these cases the formation of a condensate 
from a quantum scalar field is not a trivial process and the necessary conditions for the existence of 
such coherent behaviour must be investigated. 

In quintessence models a classical scalar field is the basic content of the model and its energy density 
is interpreted as dark energy. Although in the framework of popular particle physics models such 
as super symmetry, supergravity, and string theory many eff^orts have been concentrated on finding 
candidate scalar fields to play the role of quintessence [U] , little work has been devoted to understand 
what are the necessary conditions for a quantum scalar field to condense in a manner which satisfies the 
very special characteristics of a quintessence field. For instance, such a condensate must initially have 
a very small density, much smaller than other content of the Universe (smallness problem). Present 
observations show that dark energy behaves very similar to a cosmological constant i.e. with the 
expansion of the Universe its energy density does not change or varies very slowly. Such a behaviour 
is not trivial. In the classical quintessence models usually the potential of the model is designed 
such that a tracking solution do exist. Potentials with such property are usually non-normalizable. 
Moreover, they don't directly correspond to potentials (or kinetic terms) expected from fundamental 
theories such as supersymmetry, supergravity or string theory. Therefore one has to relate them 
ad-hocly to some sort of low energy effective model of a fundamental theory. 

The purpose of the present work is to fill the gap between quantum processes producing various species 
of particles/fields in the early Universe - presumably during and after reheating - and their classical 
component as defined in In another word, we want to see how the microscopic properties of 
matter is related to macro-physics and vis versa. We are particularly interested in condensation at 
very large scales, relevant to dark energy models. For other phenomena such as baryo- and lepto- 
genesis and Higgs mechanism, if the energy scale was much larger than Hubble constant of the epoch, 
the process can be studied locally. This is not the case for dark energy which seems to be uniform at 
largest observable scales, and therefore the expansion of the Universe could play important role in its 
evolution. 

As the quantum physics of the epoch just after reheating is not well known, we consider a simple 
toy model with a light scalar as quintessence field in interaction with two heavy scalar fields. Our 
aim is to study the evolution of the classical component - the condensate - of the quintessence field. 
Between many possible types of quantum scalar field and interaction models, we concentrate on a 
class of models in which the heaviest of three particle decays to other fields. The motivation for 
such a model is the results obtained from the study of the effect of a decaying dark matter on the 
equation of states of the Universe [7]. It has been shown that a FLRW cosmology with a decaying 
dark matter and a cosmological constant behaves similar to a cosmology with a stable dark matter 
and a dark energy with w = P/p ^ — 1. This is effectively what is concluded at least from some of 
present supernovae observations [5]. More recently, the same effect has been proved to exist for the 
general case of interaction between dark matter and dark energy [9]. It has been also shown |10] that 
if a decaying dark matter has a small branching factor to a light scalar field, the observed density 
and equation of state of dark energy can be explained without extreme fine-tuning of the potential 
or coupling constants. In other word, such a model solves both smallness and coincidence problems 
of dark energy. The present work should complete this investigation by studying the formation and 
evolution of classical component from a quantum point of view. More generally, it is believed that all 
the particles are produced directly or indirectly from the decay of inflaton or curvaton (in curvaton 
inflation models) oscillation. Quintessence field is not an exception and irrespective of the details 
of its physics, it has to be produced from the decay of the inflaton or another field. Although the 
toy model considered here basically assumes a long life heavy particle, in each step of calculation we 
also mention the differences in the results if the life time of the decaying particle is short. The main 
difference between these cases is the time duration in which the production of quintessence scalar by 
the decay is significant. 
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In Sec. [2] we describe the Lagrangian of a decaying dark matter model and evolution equation of 
the condensate. We consider three decay modes for the heavy particle and use the closed time path 
integral method to calculate the contribution of interactions in the condensate evolution. The same 
methodology has been used for studying inflation models jllj . late-time warm inflation [12], the effects 
of renormalization and initial conditions on the physics of inflation |13j , baryogenesis [18] , and coarse- 
grained formulation of decoherence |21j (see also [23] and references therein). In SeclSjwe solve field 
equations and discuss their boundary conditions. In Sec. [Hwe obtain an analytical expression for the 
asymptotic behaviour of the condensate and discuss the importance of the quantum corrections. We 
summarize the results in Sec. [SJ In Appendix |X] we obtain non- vacuum Green's functions in presence 
of a condensate. In Appendix iBl we generalize the description of a condensate to a system in which 
not all the particles are in the ground state. Appendix O presents the solution of evolution equations 
in matter dominated era. Finally, in Appendix [Pl propagators in a fluctuating background metric are 
determined. 

2 Decay in an Expanding Universe 

We consider a simple decay mode for a heavy particle X to a remnant that includes only 2 types 
of particles: a light scalar $ - light with respect to decaying particle X - and another field A of an 
arbitrary type. In fact, in a realistic particle physics model, most probably A will not be a final stable 
state and decays/fragments to other particles. Therefore it should be considered as an intermediate 
state or a collective notation for other fields. In the simplest case studied here all the particles are 
assumed to be scalar. Extension to cases where the decaying particle X and one of the remnants are 
spinors is straightforward. The quintessence field $ however, must be a scalar. We do not consider 
the condensation of vector fields here. In the extreme density of the Universe after reheating, apriori 
the formation of Cooper-pair composite scalars from fermions is also possible. This process needs 
a relatively strong interaction between fermions and can arise in local phenomena such as Higgs 
mechanism and leptogenesis which happen at high energies (short distances) [TU], but not for dark 
energy which must have a very weak interaction and acts at cosmological scales. 

We consider the following decay modes: 



Diagram ([Sj-a) is the simplest decay/interaction mode. Diagram ([Sj-b) is a prototype decay mode when 
X and share a conserved quantum number or A and A (here A = A\s considered) has a conserved 
quantum number. For instance, one of the favorite candidates for X is a sneutrino decaying to a 
much lighter scalar field (e.g. another sneutrino) carrying the same leptonic number |26j [27j . With 
seesaw mechanism in the superpartner sector (or even without it [27]) if SUSY breaking scale is lower 
than seesaw scale, a mass split between right and left neutrinos and sneutrinos will occur. As the 
right-hand neutrino super-field is assumed to be a singlet of the GUT gauge symmetry, it has only 
Yukawa-type of interaction. In such a setup X can be a heavy right sneutrino decaying to a light 
sneutrino with the same leptonic number and a pair of Higgs or Higgsino [28] . In place of assuming 
two A particles in the final state we could consider them as being different A and A, its anti-particle. 
But this adds a bit more complexity to the model and does not change its general behaviour. For 
this reason we simply consider the same field. Diagram ([Sj-c) is representative of a case where X and 
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A can be fermions (although we do not consider this case here) and <I> can be complex and carries a 
conserved charge |29j (again for simplicity we do not consider this case here either). 



The corresponding Lagrangians of these effective interactions are the followings: 

A 



int 




n 



-g^'^d.Ad.A - -m\A^ - -A"-' 

g^XA, For ©-a 
'g{g^XA^, For ©-b 
g^'^XA, For ©-b 



(6) 
(7) 
(8) 

(9) 



In addition to the interaction between X, $, and A we have assumed a power-law self-interaction 
for $ and A. If A is a collective notation for other fields in more realistic models, its self-interaction 
corresponds to the interaction between these unspecified fields. Again for the sake of simplicity in the 
rest of this work we consider A' = 0. The unstable particle X is assumed to have no self-interaction. 
Note that the same Lagrangian can be considered to present the interaction between these fields. 
Therefore with little modification, the results of this work become applicable to the case of interaction 
between dark matter and dark energy. 

Although the model presented here is quite general, we are primarily interested on the physics of 
dark energy. In this context, the heavy particle X is a candidate for the dark matter, and $ is the 
quintessence field, and the Lagrangian ([9|) presents the interaction between these fields. It is necessary 
that X and have a very weak interaction with each other and with the rest of the Universe presented 
by A. Therefore couplings A and g must be very small. 

In a realistic particle physics model, renormalization as well as non-perturbative effects can lead to 
complicated potentials for the scalar fields. An example relevant to dark energy is a pseudo-Nambu- 
Goldston boson field. Its potential is assumed to have a shift symmetry [3D]. This class of models are 
interesting for the fact that the mass of the quintessence field does not receive quantum corrections 
and can be very small. Moreover, they can be easily implemented in SUSY theories along with right- 
neutrinos and sneutrinos, as candidate for X |20j . The power-law potential considered here can be 
interpreted as the dominant or one of the terms in the polynomial expansion of the potential. In 
addition, we Will see in Seci3]that only the few lowest order in the expansion play a significant role 
in the late time behaviour of the condensate. 

The general aspects of the analysis presented here do not depend on the details of the particle. Here 
our aim is an analytical investigation of the evolution of the condensate to see whether it is possible 
at all to have a quantum condensate at large scales. To achieve this goal we had to apply many 
approximations and simplifications. A more precise solution needs numerical analysis and we leave it 
to a future work. 



2.1 Decomposition 

We decompose ^{x) to a classical (condensate) and quantum components: 

^{x) = ^{x)I + <j){x) {^) = {^^^) = ^{x) (0) = (^l^l^r) = (10) 

where / is the unit operator. Note that in (jlOp both classical and quantum components depend 
on the spacetime x. In studying inflation it is usually assumed that the pre-inflation Universe was 
homogeneous or the very fast expansion of the Universe had washed out all the inhomogeneities and 
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the condensed component of inflaton became homogeneous. Here we are studying the evolution of a 
scalar field after inflation when the distribution of the unstable X particles can have non-negligible 
inhomogeneities, specially if the decay is slow and perturbations have time to grow. Thus, we don't 
assume a homogeneous Universe, but anisotropics are assumed to be small. Moreover, for solving 
dynamic equations, in some situations we have to ignore anisotropics all together, otherwise the 
problem would not be tractable analytically. 

We assume {X) = and {A) = 0. Justification for these assumptions is the large mass and small 
coupling of X and A which should reduce their number and their quantum correlation. In other words, 
when mass is large, the minimum of the effective potential for the classical component is pushed to 
zero (see (jlSMlSp below). We find a quantitative justification for negligible condensation of massive 
fields in Sec. [H As X and A have a very weak interaction with (/> and the condensate, their evolution 
can be studied semi-classically by simply considering the decay and interaction cross-sections. Such a 
study has been already performed in [3 |10]. Therefore, here we concentrate on the evolution of the 
condensate ^{x) and if necessary, we use some of the results from the works mentioned above. 

The Lagrangian of $ is decomposed to: 



d'^x^J—g 



-a 



-9 < 



■m%ip(j) 



, n— 1 

-E 



n 



1=0 



gipXA + g(pXA For P-a 

gipXA^ + g(j)XA^ For ©-b 

^ g^'^XA + 2g^(t)XA + g(j)'^XA For ©-c 



n 



[II] 



+ 



(12) 



The first two terms in Cint are obtained after integrating out two total derivative terms. Lagrangians 
and C(j) are the same as ([6|) with <I> — )■ 99 and <1> — )• respectively. The self-interaction terms in (jl2p 
include a term proportional to (j)^ that contributes to the mass of quantum component. Therefore, 
when we use the free Lagrangian of the quantum component to compute quantum corrections, this 
term is considered to belong to C^. 

The evolution equation for the condensate (classical component) (p can be obtained from Lagrangian 
-L<j> by variation method, we must also take the expectation value of operators on the state 



1 



-My 



A 

^g^'^d^^) + + - V(i + 1) 
n ^-^ 

i=0 



^g^"'d^<p) + ml^ + - V(i + 1) 
n ^-^ 

i=0 



A 



n 



n-l 

E 

1=0 



n 
i + l 



n 
i + l 



n 
i + l 



in—i—l\ 



^'{r-'-')-g{XA)=0 

For ©-a (13) 

vHr-'-')-g{xA') = o 

For ©-b (14) 
2gip{XA) - 2g{^XA) = 

For dSD-c (15) 



Note that in ()12p non-local interactions, i.e. terms containing derivatives of ip do not contribute 
in the evolution of (p because they are all proportional to (p. After taking expectation value of the 
operators they cancel out because {(p) = by definition. The expectation values depend on the 
quantum state of the system 1^*) which presents the state all the particles in the system. From the 
structure of Lagrangians (|lip and (|12p it is clear that the mass of quantum component (p and thereby 
its evolution depends on ip. Moreover, through the interaction of $ with X and A the evolutions of 
all the constituents of this model are coupled. See Appendix [X] for more details. 
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For n ^ 2 the expectation values (c^'^""*^^)) modify the mass and self-couphng of y?. Another important 
observation is that in general the form of the potential for the classical field (p is not the same as the 
potential in the original Lagrangian, although they have the same order. Therefore, the usual practice 
in the literature of using the same potential for both quantum and classical component is only an 
approximation. 

For models (a) and (b) the expectation value of interaction in (|13p and (|14|) contains only X and A. 
But when these terms are expanded, see equations ([22]) and ([23]) below, they depend also on ip. We 
will see later that these terms play the role of a feedback between production and evolution of the 
classical field. In particular, they prevent a complete decay of the condensate with the accelerating 
expansion of the universe. This effect is similar to what was found in the classical treatment of the 
same models in Ref. [10]. Interaction (c) is more complex and various evolution histories for ip are 
possible. They depend on the value and sign of g the coupling of <I> to X, self-coupling A, and the 
order of self-interaction potential n. For instance, the mass can become imaginary (tachyonic) even 
without self-interaction and lead to a symmetry breaking. Tachyonic scalar fields have been suggested 
as quintessence field specially in the framework of models with w < —1 [31]. We discus the difference 
between these decay modes at each step of calculation. 



2.2 Expectation values 

We use Schwinger closed time path formalism also called in-in to calculate expectation values. Recent 
reviews of this formalism are available [22] [23] and here we only present the results. Zero-order (tree) 
diagrams for the expectation values in (jl3m5|) are shown in (|17p . (|18p and p9p . The next relevant 
diagrams are of order and for dark energy models are expected to be negligibly small. Evidently 
the decomposition of $ also affects the renormalization of the model. This issue has been already 
studied [13] and we do not consider it here. One example of higher order diagrams is shown in ()16p . 
These types of graphs are specially important for studying renormalization in the context of a realistic 
particle physics model. Thus, for the phenomenological models considered here we ignore them. 



{XA)a = 
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The index a and c in ()16p and (jlSp refer to the corresponding interaction model respectively. The 
graph ()20p is an example of self-interaction terms in (jl3H15p with n external lines where i ^ 1 of them 
are of type (p and the rest of type ip. The dash lines present the classical component (p. When i = 1 
there is an additional interaction involving derivative of the classical field. 



(0) 



-0 



+ 



-• 

s 

\ 



(21) 



At lowest order the sum of these graphs is null because they correspond to the dynamic equation of 
see (|13m5p . This is consistent with the decomposition ([TO]) . 



The corresponding expectation values at zero order are: 



{XA)a 

{XA^) 
{(pXA) 



-iX 



^d'^yfiy) G^{x,y)G^{x,y)G^{x,y) - G^{x,y)G^{x,y)G'^{x,y) 
^dWiy) G>{x,y)G>ix,y)-G<ix,y)G<{x,y) 
~gd^y^{y) G>[x,y)G>{x,y)G>{x,y)-G<[x,y)G<{x,y)G<{x,y) 
-gd'y^--\y)\[G>{x,y)r -[G<{x,y)Y 



Advanced and retarded propagators G^ and are defined as: 

G^{x,y) = -i{'>p{x)'>p\y)) = -itr{ip{x)Tp''{y)p) 
G^{x,y) = ^i{i^\y)i){x)) =^iii{il)\y)il:{x)p) 



(22) 
(23) 
(24) 
(25) 
(26) 



(27) 
(28) 



where ip{x) presents one of (j), X ox A fields and p = |^')(^| is the density (projection) operator for 
the state \^). The upper and lower signs in (j28p are respectively for bosons and fermions. Definitions 
(|27p and (|28p correspond to the general case of a complex field. Here we only consider real fields and 
therefore ip{x) = iIj'^{x). Feynman propagators are related to G^{x,y) and G^(x,y): 

GF{x,y) = -i{Ti;ix)ijHy)) = G>{x,y)eix^-y') + G<{x,y)e{y'-x') (29) 
GF{x,y) ^ -i(TV(x)Vt(y)) = G>(x,y)e(2/°-xO) + G<(x,y)e(xO-2/0) (30) 

The next step is the calculation of propagators. 



2.3 Propagators and the evolution of quantum components 

Feynman propagators Gp{x,y), i = (p, X, A can be determined using field equations from La- 
grangians ([6])-([9]). The free equations of motion lead to the following equations for the propagators: 



1 



-A^iV^g^^d^G^x - y)) + {mi + (n - l)Av."-2)G^(x - y) 



-AiV^a^'d^GU^ - y)) + m^GUx - y) 



.5^{x - y) 



.5'^{x - y) 



X.A 



(31) 
(32) 



The free propagator of (j) is independent of the type of interaction between X, A, and and therefore 
equation (j3T]) is valid for all the interaction models in ([5]). Note that G'p{x — y) is coupled to the 



8 



condensate field ip even at classical level. On the other hand, evolution equations ()13p -()15 p depend 
on the interaction between quantum fields X and A. This means that all the components of 
the model are coupled even at lowest order. The coupling between quantum component (j) and the 
classical component (/? is the origin of the back-reaction of the condensate formation on the quantum 
fields. Its strength depends on the mass, order, and strength of the self-interaction. Assuming that 
initially $ particles are produced only through the decay of X, the initial value of (/? = 0. Therefore 
initially the coupling between (f) and ip was very small. With the growth of the (p amplitude, the 
effective mass of (p particles becomes larger than their bare mass. On the other hand, this affects the 
growth of the condensate because due to an energy barrier (p particles will not be able to join the 
condensate anymore. Therefore, there is a negative feedback that prevents an explosive formation of 
the condensate. 

Assuming a quick, roughly immediate decoherence of (f) and other speciesH, if their interaction is weak, 
we do not need to consider a complete quantum treatment of this modeo, and therefore the quantum 
state of free particles |^'/), including (j), can be approximated by direct multiplication of single particle 

l^/) = E(2)/^(^'{^'^»l?'P (33) 

The indices i and j present the species type and particle number respectively, and {pj} the momentum 
of all states. In the context of dark energy models single particle description is a good approximation 
because the interaction of ^ with itself and with other particles must be very small. This can be not 
generalized to process in which fields can have strong couplings, such as supersymmetric models before 
SUSY breaking and electroweak in the early Universe, QCD, etc. In this case a complete N-particle 
evolution of the unfactorizable wave functions must be considered [U [23] • 

In the Introduction we argued that the expectation value of <I> on the state |^/) is null, f\^\^ j) = 0. 
The complete state of the Universe l^*) depends on both condensate and free particles i.e. = {"^f, ip). 
When the coupling between these components is small ^ can be factorized: 

\^) = \^f,ip) f^l^f) (S)\ip) (34) 

A special description for | ip) is given in Appendix [A] and a more general one in Appendix [Bj The 
presence of in ^'j in the definition of 1^*) reflects the fact that because of their interaction we 
can not completely separate these subsystems. The amplitudes of single-particle states /* - the one- 
particle distribution functions - in ([55|) are considered to depend on the spacetime coordinates to 
reflect the process of squeezing and decoherence of the wave functions. In fact in this setup, the 
only difference between squeezed (classical) particles of the same species is their place with respect 
to the random fluctuations of the background. More precisely quantization of inflaton and other 
fields induces a quantized metric fluctuation - corresponding to a semi-classical treatment of gravity - 
because the metric is related to matter through Einstein equations (see e.g. [5] and references therein). 
The decoherence of inflaton oscillations and other fields which are produced by the decay of the 
inflation oscillations make particles to behave classically. At the same time this process decoheres the 
metric fluctuations because metric is treated as a secondary (dependent) field. The result is a classical 
distribution of free particles but the particles/fields which stay correlated. Therefore we can interpret 
/* as a classical distribution. 

In the classical limit the evolution of /*(x,p) is governed by the Boltzmann equation |24j : 

8 f (*) 

P^9^/«(^,p) - (ri;,pV)^ = -iA{x,p)+B{x,p,^))f(-'Hx,p) + cix,p,^) + 

V{x,p,ip) +£{x,p). (35) 

'^This assumption seems inconsistent. In one hand we want tliat $ has a very weak interaction with itself and with 
other particles. On the other hand we simplify the problem by assuming that they decohere quickly. We must also 
remind that if X particles that produce $ are decohered, so do $ particles. 

^For a complete treatment one has to use the techniques of non-equilibrium quantum field theory and Kadanoff-Baym 
equations (see for instance [4] [23]) which lead to quantum Boltzmann equations (see [17] and references therein.) 
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A{x,p) = Tirrii. (36) 

C{x,p) = T.^jrn,j^J dp,f'-^\x,p,)^^. (38) 

^(^'^') = Et^^ / dp,dpkf^^\x,p,)f^PHx,pk)A{s) '^'''+p^- . (39) 

MPuPj) = {{Pi-Pj? - mjm])^ = -mf- m]f - Amjm])^. (40) 



where m^, Fj = I/tj, Tj are respectively mass, decay width, and hfetime of species i; cjjj is the total 
interaction cross-section between species i and j at a fixed center of mass energy s; daj+k^i^„Jdp = 
Eda / Qip'^ dpdO, is the Lorentz invariant differential cross-section of production of i in the interac- 
tion of j and k; gi is the number of internal degrees of freedom (e.g. spin, color, etc.); dp = d^p/E is 
the Lorentz invariant measure in momentum space; the term dM^^^ j/dp is the differential multiplicity 
of species i in the decay of j; and finally T^p is the Levi-Civita connection. Note that the right hand 
side of (|35p is written in the local Minkovski frame in which the expression for the cross-sections 
is simple. As the effective mass of (f) depends on the classical field the right hand side of (j35p 
and thereby distributions P{x,p) depend on ip. Thus, as mentioned at the beginning of this section, 
evolution of quantum and classical components of this model are coupled. 

Finally we must add Einstein equations to the set of evolution equations discussed above. We only 
consider scalar fluctuations in the linear regime and assume that the deviation of Z*-*^ from a perfect 
fluid is small. With these simpliflcations the metric in Newtonian gauge can be written as: 

ds^ = {l+2i;{x))dt'^-a'^{t){l-2'ilj{x)5ijdx'dx^ = a'^{7])[{l+2'ilj{x))dif -{l-2ip{x)5ijdx'dx^)], dt = adi] 

(41) 

where t and r] are respectively comoving and conformal times. Einstein equations for this metric can 
be found in textbooks, see for instance [3]. We do not reproduce them here because in the present 
work we do not solve them along with held equations. Nonetheless we only remind the evolution 
equation for a{t) which its evolution has special importance for dark energy: 

ct^ 8ttG v-^ , , , , , 

i 

where pi is the energy density of species i. During radiation domination epoch the density of non- 
relativistic particles such as X is by definition negligible, and evolution of a{t) is governed by relativistic 
species which are not considered here explicitly. From the observed density of dark energy we can 
conclude that in this epoch its density was much smaller than other components, and had negligible 
effect on the evolution of expansion factor. In the matter domination epoch both X and A are assumed 
to be non-relativistic. If the lifetime of X is much shorter than the age of the Universe at the beginning 
of matter domination epoch, most of X particles have decayed, and it does not play a significant role in 
the evolution of a{t) which is determined by other non-relativistic species. If the lifetime of X is much 
larger than the age of the Universe, then X particles can have a significant contribution in the total 
density of matter. As it decays very slowly, in the calculation of a{t) it can be approximately treated 
as stable. In this case the evolution of a{t) would be similar to a CDM model. A better estimation 
of a{t) can be obtained by taking into account the decay of X to relativistic particles. This method 
is used Ref. [7j and a(t) is calculated. However, here we use the simpler approximation because the 
problem in hand is very complex and we want to keep a{t) decoupled from other equations. At late 
times when the density of condensate becomes comparable to matter density the full theory, including 
Boltzmann equations must be solved. In this case the evolution of a(t) is not simple and needs a full 
numerical solution and we leave it to a future work. 
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Quantum interactions happen at high energy scales, i.e. short distances. Thus, in the dynamic 
equations of the quantum fields we use only the homogeneous metric. This approximation is valid if 
at the epoch just after the production of X particles H <^ max- Assuming a low energy preheating 
temperature ~ 10^ GeV and mx ~ ^gut ^ 10^^ GeV, mx would fulfill the above condition. However, 
at early times this condition is not satisfied by the self-interacting field ^> because its mass is expected 
to be very small. In Appendix ([D]) we show that at linear order the propagator of the fields in a 
background with small fluctuations are simply Gfi{x, y){l + ^l^). Therefore we can use the homogeneous 
background to solve dynamic equations and then correct it for the effect of small fluctuations if 
necessary. For this reason in the following we mainly use a homogeneous background and discuss 
the effect of fluctuations afterward. Apriori the evolution equation of the condensate also must be 
written for a fluctuating metric, specially if $ presents the quintessence field of a dark energy model. 
In Appendix ([D|) this equation is determined in Newton gauge. However, due to the non-linearity of 
this equation, it is not possible to find a simple perturbative correction of the solution. Therefore in 
the present work we neglect metric fluctuations in this equation. 



3 Solution of the evolution equations 

In the previous section after writing all the dynamic equations it became clear that they are coupled and 
in general non-linear. In this situation it is impossible to proceed analytically and obtain a solution 
for these equations unless we break the mutual coupling of equations by taking some simplifying 
approximations. This should be possible because both mass and couplings of <I> are considered to 
be small. Therefore, the contribution of the condensate ip in the evolution of particle distributions, 
equation ()35p . must be small. Moreover, in this work our focus is on the formation of the condensate 
and an approximate distribution for other particles should be enough for a zero order estimation of 
the condensate evolution. Therefore, in place of solving the complete set of Boltzmann equations, we 
consider initial thermal distributions for the particles and assume that the interaction terms in the 
right hand side of leads to a slight difference between the effective temperature of species. This 
simplification is similar to temperature shift considered for the CMB in the treatment of SZ effect. 
Giving the fact that uncertainties on the initial temperature of the species is large, a slight modification 
should not have large effect on the condensation process and the properties of the condensate, otherwise 
the model would need fine-tuning and lose its reliability. 

The advanced and retarded propagators defined in (j27p and ()28p . are needed for the calculation of 
expectation values in (j23ti25p . They are calculated on the non- vacuum states of the matter content of 
the Universe. Appendix Rl explains in details the derivation of such propagators from vacuum Green's 
functions. Energy-momentum distributions of particles discussed above are necessary at this step. 
They make the connection between macroscopic cosmological evolution of the distribution of species 
and micro-physics of the condensate. 

To solve the Green's functions and the evolution equation of the condensate, it is more convenient to 
write them with respect to conformal time rj along with a redefinition of Lp: 

X = a^ (43) 

Then, if we considering only the homogeneous component of metric ()4ip . the evolution equation of 
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the classical field x for the three decay models in ([5]) are the followings: 



X" - S^AdjX + {a'rnl - + - ^ a^-\i + 1) . J xHr^'~') - ah{XA)a = 

For (I5])-a(44) 

X" - ^.AdjX + {a'ml - + - ^ a^-\i + 1) . J X^(<^""*"') - a3g(XA2) = 

CL Til \ Z ~i i / 



For (l5])-b(45) 



x!' - kAQ^X + - + - J] a3-(i + 1) . "^J x\<^"~''~^) - Ho'xkXA), - 2ga\^XA) = 

For (ISI)-c(46) 



Note that even if we neglect the expectation values related to the interaction between X, A and due 
to the interaction of classical component (the condensate) with the quantum component the effective 
potential of the condensate (/? is not equal to the self-interaction term in the original Lagrangian and 
includes quantum corrections. In Sec. H] we discuss the circumstances in which these corrections play 
an important role in the evolution of ^p. 

The vacuum propagators of quantum fields 0, X and A areH: 



^Gl[x, y) - 6ijdid,Gj{x, y) + {a^ml -^ + {n- l)Aa^-"x"~')G'J(x, y) = -i *^'^^^ ^1 ^7) 

^GUx, y) - 6ijdidjGUx, y) + {a^nl - y) = i = x, A (48) 

Gl = a{r^)G'^p , G\ = a{r^)G^ , G'p = a{ri)Gi (49) 



with /' = df/dr]. Note that A, X and T indices are respectively another name for A, X and (f) used 
only in the modified propagators which are defined in (|49p . The classical field x (or equivalently f) 
contributes in the mass term of the quantum component T (or equivalently (p) . This term couples the 
Green's function equation of (p to the evolution equation of condensate i.e. to one of the equations 
()44M6p depending on the decay model. As the coupling A is assumed to be small, we can linearize ([17|) 
and use a WKB-like prescription to obtain an approximate solution. However, for solving evolution 
equations ()44M6p the linearization is not always a good approximation. In Sec. [5] we will argue that a 
late time non-zero slowly varying condensate can be obtained only when the full non-linear equations 
are considered. 



3.1 Homogeneous solution of field equations 

The X particles are presumably produced during reheating epoch |25j and their decay begins afterward. 
In this epoch relativistic particles dominate the energy density of the Universe, thus we first consider 
this epoch. Fortunately, in this epoch the homogeneous field equations have exact solutions. In 
matter domination epoch only for special cases an analytical solution exists. They are discussed in 
the Appendix O 

The expansion factor a{7]) in the radiation dominated epoch has the following time dependence: 

a = ao( — ) = qq— a" = (50) 

^Moat of the calculations in this work are performed on a space-like 3-surface. However, to simplify the notation, we 
omit the vector sign on 3- vectors except in places where this can make a confusion. 



12 



After taking the Fourier transform of the spatial coordinates and neglecting the (/^-dependent term in 
([I7|) , the solutions of the associated homogeneous equation of (|17|) and (jlS]) are well known [32] [33] : 

{-^ + k' + a'm])Ui{v) = (51) 

i^k = J d^xU\xy^-^ i = ^ , X, A (52) 
= ciUkiiai, z'e'l) + dlVk{iai, z'e'i) (53) 

z = Oi — , 6i = V2oo%"ii = \-Fr- = - = (54) 
r/o V -no ZaQTrii Znii 

where ao and Hq are respectively the expansion factor and Hubble constant at initial conformal time 
r]Q. The functions Uk and are two independent parabolic cylindrical functions [3?j : 

C/(a,z) = 2/icos7r(^ + |) -y2sin7r(^ + ^) (55) 
V{a,z) = yisin7r(^ + ^) + y2Cos7r(^ + ^) (56) 

yi{a,z) = i 1 ' e 4 ,Fi(- + -;-;-) (57) 



4 ' 2 



^4 2' 2' 2 



V^2 

y2(a,.) = %^e-x,Fi(^ + ^;^;^) (58) 



-/7r24 ' 2 



4 2' 2' 2 



Prom now on for simplicity we drop the species index i except when its presence is necessary. We 
call two independent solutions of ([5T]) in a general basis Uk and Vk- If we want that these solutions 
correspond to the coefficients of the canonical decomposition of (p, equation (jl02p in Appendix [Aj we 
must choose a basis such that Vk =^k- ^^^^ ^^^^ ^^^^ work we only consider this basis. Note that 
the two solutions U and V in (|53p are not complex conjugate of each other and therefore can not be 
identified with and Vk- 

We are interested on the asymptotic behaviour of U and V functions when rj/rjQ ^ 1. Their asymptotic 
expressions are (see e.g. [35] and references therein): 

C/.(za.,.Vf) = e^-^i.'-^e-^(-'■^^'+V)^(2^)^(i + !|i), ^±i2i (59) 

4 2 slz' s 



Vk{iai,z'e'^) = iUk{iai, z' + \ -T{-iai + -)e 4 'sz' 

V vr z 

E(-20-(^). (60) 
The evolution equation for free Feynman propagators - the 2-point Green's functions - is defined as: 



+ e + a'm^)Gk{rj, rj') = (62) 



When 7] ^ rj' , equation (j62p is the same as the homogeneous equation ()5ip . and therefore the solution 
of ()62p is a linear combination of two independent solutions of ()5ip . According to the definition 
of Feynman propagators (j27p and (j28p . they can be divided to advanced and retarded propagating 
components and . The transformation r] ^ -q' changes the role of these propagators, -f-)- G^. 
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After adding the effect of a non- vacuum state, as explained in AppendixjAj the non- vacuum propagator 
G{r], r]') has the following expansion: 



Q{rj - 7]') + 
6(7/' - r/) 



(63) 



where , 



and B< 



are integration constants. In the Appendix |X] we show that for the 
free propagators - at the lowest perturbation order - if the state |^') is not vacuum, it is possible to 
include its effect in the boundary conditions imposed on the propagator. Comparing (j63p with ()106p 
in Appendix Rl the relation between these constants and the initial state can be found: 



kik2...kn 



1^' 



k\k2...kn I 



(64) 
(65) 



See Appendix lAl for a detailed description of It is easy to see that with relations and ([65 
between constant coefficient, the consistency condition defined as: 



ri=ri' 



(66) 



is automatically satisfied. Therefore propagators over a non-vacuum state depend only on this state 
and the solutions of the field equation. They depend also on two arbitrary constants c\ and in the 
solution of field equation. They must be fixed by the initial conditions too. 

There is one more consistency condition that propagators must satisfy. By integrating the two sides 
of the equation ()62p with respect to r] in an infinitesimal region around 77' and by using the solution 
([63|) . we find the following constraint: 



(67) 



This relation fixes one of the integration constants in (j53p . It is easy to see that a constant shift of 
the argument of Uk (??) does not violate ()67p . This means that the phase of Uk (rj) is not an observable 
and can be fixed arbitrarily. It is also interesting to note that multiplication of two sides of (j67p 
with an arbitrary constant rescales a{rj) which is equivalent to redefinition of oq. Rescaling of cq is 
equivalent to redefinition of coordinates and therefore is not an observable. In Minkovski spacetime 
the scale factor a is fixed to 1, and therefore there is no place for rescaling. In another word, in a 
Minkovski space the normalization of the propagators is an observable and affects the final results. 
This scaling property in FLRW and De-Sitter metrics is a consequence of diffeomorphism invariance 
in the framework of curved spacetimes and general relativity. 



3.2 Initial conditions for propagators 

Field equations are second order differential equations and need the initial value of the field and its 
derivative or a combination of them to obtain a complete description of the solutions. The general 
initial conditions for a bounded system - including both Neumann and Dirichlet conditions as special 
cases - are the followings [38] : 

n^'df.U = -iKU , g^.^n^'n" = 1 (68) 

The 4-vector is the normal to the boundary surface. If the boundary is space-like, can be 
normalized as = a^^(l, 0, 0, 0). Then: 

a-^d^U = -iUA (69) 
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The constant fC depends on the scale k. In any boundary problem, the boundary conditions must 
be defined for all the boundaries. Thus, in a cosmological setup the initial condition constraints (j69p 
must be applied to both past (initial) and future (final) boundary surfaces [38] • But in the case of 
propagators, they are respectively applicable to past and future propagators only. In each case the 
other boundary condition is the consistency condition ()66p . Assuming different values for /C on these 
boundaries, we find: 

^^• = ^4t^ ' ^■ = ^'^ (70) 
ajUkivj) 

Indexes i and / refer to the value of quantities at initial and final 3-surfaces. These boundary conditions 
relate /Cj and fCf to and dk in ([55]) . In fact using ([U7|) along with we find: 

i^^(^^-^'' = a2(r?,)(/c,(fc,^-)+/c*(fc,^,)) ' '^^(^^■)'' = ^^ifcl^yrl^^ ' = ''^ (^^^ 

Application of (j7ip to the solution ()53p at two boundaries fixes dynamical constant and as 
functions of /Cj and /Cj up to a constant phase. In fact as we mentioned previously, due to the equality 
(j67p the phase of is not an observable. Therefore, we can assume that it is zero on the initial and 
final boundaries. In the next section we will calculate propagators and evolution equation of x 
radiation and matter domination epoch separately. Thus, the initial and final boundaries correspond 
to the beginning and end of each epoch, and we must respect continuity condition, i.e. the initial 
condition for one epoch corresponds to the end condition for the previous epoch. 

In a cosmological context ICf can be decided based on observations, but /Cj is unknown and leaves 
one model dependent constant that should be fixed by the physics of early universe and the special 
state of our Universe among all possible states. This arbitrariness of the general solution or in other 
words the vacuum of the theory is well known [iQ]. In the case of inflation, a class of possible vacuum 
solutions called a-vacuum are usually used: 

ICi,JCf = ^k^/alj + m^ (72) 

and one obtains the well known Bunch-Davies solutions 1381. 



Alternatively one can fix the solution Uk at one of the boundaries and apply the boundary condition 
(|70p only to that 3-surface. Although this does not solve the problem of arbitrariness of /C and its 
k dependence, it reduces it to only one of the boundary surfaces, for instance to the final 3-surface. 
This make the choice of (|72p physically motivated. In addition, by applying (|7Up only to one of 
the boundaries, the causality of the constraint ()69p is more transparent and the state of the second 
boundary is directly related to the physics of the first one through the evolution equation. We remind 
that the evolution of the expansion factor 0(77) is related to all type of matter including the condensate 
through the equation (jl2]) . 

An additional and somehow hidden arbitrariness in this formalism is the fact that apriori k dependence 
of the boundary constant /C does not need to be the same for all the fields of the model. However, 
different k dependence breaks the Equivalence Principal. Similarly, a value different from (j72p for fC 
will lead to the breaking of the translation symmetry [38] [13] . In the context of quantum gravity the 
violation of both these laws are expected and therefore, in a general framework the choice of different 
fC for fields is allowed. 

Before finishing this section we discuss also /c-dependence oi and T4. This is specially important 
in the context of dark energy because no strong fluctuation is observed in its spatial distribution. 
The solutions of equation ([5T|) depends on k only through defined in From the asymptotic 

expression oiUk it is evident that apart from a constant term which can be absorbed in the integration 
constants, other terms containing a, are either oscillating or approach a constant for ^ 1 which 
is equivalent to |A;| — )• 00. It can be shown that Vk also has the same type of behaviour because 
|r(— ioi + l/2p = vr/ coshTTOj. As for the integration constants and d^, from ([7T|) and ([72]) it is easy 
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to see that for large \k\ they are proportional to Thus the amplitude of the Green's functions 

is asymptotically proportional to l/\k\ and an oscillatory k dependent component. 

3.3 Propagator of light scalar 

Finally after finding the solution of free propagators ()47M9p without f dependent terms, we use a 
WKB-like approximation to correct the propagator of (p for the contribution of self-interaction in its 
mass. We neglect the effect of varying mass on and replace z' in U with: 

z' ^ 72^^ / d{^) f 1 + (n - 1)A^^^') ' (73) 

where ^pk is the Fourier transform of (p{x) (see Appendix lAl for the details of approximations) and 0$ 
is defined in (j54|) . 

With this correction we have the solution of all the free Green's functions at lowest order. However, 
the dependence of on (/j(x) (or equivalently x) couples propagators to the evolution equation of ip 
i.e one of the equations (j44p . (j45p or (|46p . depending on the interaction mode. A complete solution 
can be obtained only through numerical calculations. Nonetheless, when we calculate expectation 
values, the correction (j73p produces high order of 99 in a polynomial expansion which at lowest order 
of approximation can be neglected. This simplifies analytical tracking of the evolution of condensate. 



4 Evolution of the condensate 

Finally in this section we study the evolution of the condensate during the cosmic time from the end 
of massive production of X particles - presumably the end of reheating - until matter domination 
epoch. Our aim is to find analytical approximations for the solutions of equations ()44ll46p which rule 
the evolution of the condensate for the three decay/inter action models considered in this work. We 
describe the solution for the model (a) which is the simplest one in details. For other two models we 
briefly explain their differences and deviations from model (a). 

We begin with the radiation domination epoch and use the results obtained in Sec. [3l Then, we con- 
sider the matter domination epoch for which the solutions of field equation are discussed in Appendix 
O Before going to details we should make a remark about the necessity of considering quantum 
corrections in the evolution of the condensate. The reason is the complexity of equations ()44ll46p and 
thereby the solutions obtained in this section. They raise a question about the necessity of considering 
such complex model and the significance of quantum corrections in the evolution of condensate. After 
all, in other contexts such as inflation, a classical homogeneous held is considered and only fluctuations 
around it are quantized and studied. We should remind that here we are considering a situation in 
which initially the scalar fleld is uniformly null. Such an initial point is the minimum of a n > 
potential which has a clear physical interpretation in many-body quantum systems and perturbative 
fleld theories. Spinorial models in which the potential has a non-zero maximum at origin induce an 
imaginary (tachyonic) mass term and can not correspond to a fundamental field and must be con- 
sidered as an effective field. The same argument applies to inverse power-law potentials which are 
present in the majority of quintessence models. In fact, in all these cases the potential is the effective 
potential of the scalar condensate not the fundamental quantum field itself - the best example of such 
models is Higgs. Therefore, when we solve the dynamic equation with this minimum point as initial 
condition, no condensate can form unless an additional (quantum) interaction is present. In the case 
of inflation, this pre-inflationary step is usually overlooked and its presence is simply assumed as an 
initial condition. Therefore, the issue of the formation of a cosmic scalar - a condensate - discussed 
here is relevant for inflation models too. 
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4.1 Radiation domination epoch 



Model (a) is the simplest case between models considered in this work. We first neglect the self- 
interaction term in (j44|) and take the Fourier transform of the left hand side of this equation (See 
Appendix lAl for technical details and approximations that have been taken): 



X" + {k^ + a^mDx + 



(27r)6 



(fkid^k2S^^\k-ki-k2) 







(74) 



This equation is linear and we can use WKB-like methods to find an approximative solution. If we 
neglect interaction terms, equation ()74p has an exact solution of the form ()53p . Because it is assumed 
that the coupling g is smalo, a WKB-like methods allows to find an approximate analytical solution 
when the interaction terms are taken into account: 



xtH^) = crC/fc(^a<i.,Wr(r?)e<i.e*t) + 4"Vfc(ia$,Wr(^)^*e^* 



Ml 



drj{ 1 + 



-9 



^g^ 



(27r)3m|a(7?) 



d^kid?k26^^\k-ki-k2) 



dr]'y 



Gi": iv, V')GI> iv, V') - Gif iv, v')Gl< (r?, r?') 



(75) 



(76) 



where the index (a) refers to the interaction model. Considering the complexity of expressions (|119p 
and (jl20p for propagators and multiple integrals in ([76]) . it is evident that the expression (j75]) is very 
involved. But, we are essentially interested in the asymptotic behaviour of Xk to see whether it grows 
with time. For this reason, in place of presenting the full expression of the integrals, we concentrate 
on growing terms, their order, and the conditions for their existence. 

Using equations ([^^ and ([UT]) . the asymptotic expression of the functions Uk and respectively, we 
find that the asymptotic expression of Xk includes terms of the following form: 



iz' + hkiz')y 



where z' , 0$, and 9^ are defined in ([MD and hk{z') = wi''\r])9^ 
large scale modes i.e. k <C for which the parameter is expected to be small. The power-law 
factor in ()77p decays as ~ l/y^- On the other hand, terms with in which ±i{z' + hk{z'))'^ , the exponent 
of the exponential factor, has a positive real part will grow exponentially and ones with opposite sign 
will decline exponentially. Thus, in radiation domination epoch the production of $ particles by the 
decay of X alone is enough for the formation of a condensate. However, due to the smallness of the 
interaction term (|76p when r age of the Universe in that epoch, the growth of the condensate can 
be very slow. 

To complete this argument we must also show that the function fkiz'), in particular its imaginary part 
is not zero. Therefore, we must determine the multiple integrals in ()76p . Considering the complexity of 
propagators ()119p and ()120p . it is evident that their calculation is long and tedious but straightforward. 
Thus, in place of presenting the complete expression of fkiz'), we only explain the general form of 
the terms it contains and their behaviour for rj/rjO ^ 1 which is relevant to the asymptotic evolution 
of condensate if. We also use the asymptotic expressions of the functions and to expand and 
calculate the propagators in the closed time path integrals. In addition, as we discussed in Appendix 
1X1 we consider that barred coordinates in the propagators (|119p and ()120p are independent from x 
and y. They present the average space-time coordinates in which these processes occur. With these 



z' < z'. 



(77) 



We are interested in the 



*In model (a) the coupling g has a mass dimension of 1. Therefore the statement about the weakness of interaction 
in l|74|) or in other word the smallness of g refers to its comparison with (fc^ + a^m|) term. 
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simplifications, thermal distributions in the propagators do not contribute in the integration over time 
in ()76p . The integral over r( includes terms of the following form: 



e*— = -(-)2+- r(2 + -,-^x2), x = ^ (78) 

where parameters a = (±a^ ± ax) and /3 = (±^^ ± ^^) (all combinations of signs are present). 
The integral in (j78p includes the contribution of both advanced and retarded propagators that cover 
distinct time domains - for a given r/, in advanced term r]' < rj and in retarded term rj' > rj. 

The next step is taking the indefinite integral over rj which includes terms of the following form: 

- r(2 + -|.^) . -- (f )-- (^)- e -0 (79) 

For this integration we have used the asymptotic expression of incomplete Gamma functions r(a, x) ~ 
^a-ig-x^ Note that the x~^ factor in the integrand is due to the a~^(ry) factor in (??). It has an 
important role in the evolution of the condensate and its presence depends on the interaction model. 
Finally, we find the following expression for W^."^(?7) when 77/% ^ 1: 

^t\'n)^V + I d^hd'k26'^'Hk-h-k2)y^A^^{h,k2,x) e^p(i{2aln^-^)] (80) 

(27r)'^m| J ^ y Arj^ J 

The factors Aap depend on a a, ctx, Oa-, Qx, and thereby on ki and k2. The parameter a also depends 
on the momentums, but it appears only in the logarithmic term which increases much slower than 
the quadratic term. Thus, its effect can be important only at large momentums. More importantly 
Af^p factors depend on the distributions of X and A particles. Each of them includes one of the 
following factors: f^^^f^\ f'-^Hl + f^^), {l + f^^^)f^^\ or (1 + /(^))(1 + /W). This is the reason 
for the inclusion of x in ()80p . As we mentioned before, we first consider x as an independent variable. 
Then, after the calculation of closed time path integrals, we identify it with ij (neglecting its spacial 
dependence). Note that due to the constant term in some of the above factors, the interaction term 
is not zero when distributions are null, i.e. when expectation values are calculated for vacuum. This 
apparent inconsistency can be solved if the initial growth rate of x assumed to be null, which is 
consistent with the concept of vacuum as the absence of any particle, see Sec. 14.31 for details. Before 
proceeding with the final expression for Uk and Vfc, we discuss the distribution of A and X particles. 

The X particles are expected to be heavy with a mass of the order of GUT scale ~ lO'^^ GeV. On the 
other hand, reheating temperature must be < 10^ GeV to prevent overproduction of gravitinos |14j . 
Therefore, we expect that at production and thermalization epoch - if X particles have ever been in 
thermal equilibrium with the rest of the Universe [151 125] - their temperature or more generally their 
kinetic energy was much smaller than their mass. Thus, we can safely consider their temperature to 
be zero. This simplifies the expression of their propagators. In this case their density evolves as: 

Pxix) = px{xo)4e~'-^, f<^Hx,p) ^ ^-^^^^6{\p\) (81) 
a'^ mx 

Using Pxix) and the effective temperature of A particles, equation piip in Appendix[Xl we find that 

at t ^ Iq: 

^ TT^MAPx{x)to I t_-to 

^bTa-{ )4e 4. (82) 

This means that if (t — to) ^ ^ during radiation domination, after a rapid rise of Ta from zero, 
see (|113p the effective temperature of A particles quickly approaches to a constant value. Therefore, 



their distribution during this epoch approaches f^^^ ~ exp{— P a\J k'^ / ol^ + "^a) '^here Pa is roughly 
constant. Due to the presence of factor the effect of scale dependent term in /(^) deer eases with 
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time. Note that the effect of gravitational growth is included only in the effective temperature ()82p 
because we didn't solve the Boltzmann equation for A. This should be an enough good approximation 
for the needs of the present work as we are mostly interested in large scale behaviour of the fields and 
distributions of particles. 

A complete analytical expression of W^"^ needs also the integration over momentums ki and k2 ■ One 
of these triple-integrals is canceled by the delta function presenting the conservation of momentum: 
k2 = k — hi. The other integral is reduced to a double integral because without loss of generality 
the z-axis can be selected to be orthogonal to the plane k — ki. Nonetheless, these integrals are very 
complex except for the terms containing f^^^. The presence of a delta functions in these terms, see 
equation (f8T]l . reduces the second triple- integrals too. Fortunately these terms are much larger than 
terms due to the contribution of vacuum which are subdominant and can be neglected. However, 
equation ()8ip shows also that f^^^ decreases with the expansion of the Universe as well as with the 
decay of X particles. In this case, at late times the vacuum terms can become dominant. Nonetheless, 
if a significant fraction of X particles persist until the end of the radiation domination epoch, the 
terms proportional to f^^'^ continue to be dominant, and therefore Wfc has the following approximate 
expression: 

ig^??^ . . f.,^ . V Prf , 
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(83) 



In Sec. 



we showed that Green's functions are asymptotically proportional to l/|/c|. Therefore, we 
'''^^ oc l/|A;p, and large \k\ modes decay quickly. This is consistent with the behaviour 



expect that 
of dark energy. 

After applying ([83|) to and in 
the following expressions: 
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The presence of a real exponential term in both independent solutions of the evolution equation and 
the phase difference between them means that in the radiation domination epoch there is always a 
growing term that assures the accumulation of the condensate, although due to the smallness of the 
coefficients which is proportional to g^, its growth can be very slow. Therefore, we conclude that 
in this regime the production of $ particles by the slow decay of X particles according to the model 
(a) is enough to produce a growing condensate. 

It is useful to compare this result with exponential particle production during preheating. In fact 
equation (j74|) has the same structure as the linearized equation for the quantum fluctuations around 
the minimum of the inflaton potential. But here the effective potential is more complex than many 
of inflationary models. In addition, it contains space-time dependent coefficients. Nonetheless, the 
general aspects of the asymptotic behaviour of these models are similar, compare Figure [T] with e.g. 
figure 3 in Ref. [16]. The reason for this similarity is the fact that exact solutions of the non-perturbed 
equation in both cases are the same, and in our case quantum corrections include Green's functions 
which are again the solution of the same type of equation with different boundary conditions. They 
are combined in a sophisticated manner - through integration which in the approximate solutions ([8] 
and ([85]) acts like a linear operation. 
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log ii/ii „ log ii/Ti „ log ti/ii „ 
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log ril-q „ log t\Iti „ log tiIt\ „ 

Figure 1: Real, imaginary, and absolute value of one of the terms respectively in Uk and Vk- We used 
a = and /3 = 100. The general aspects of these functions are not very sensitive to a and are very 
similar for large (3 > 10. Note that although there are resonant jumps in the value of U and V, due 
to the complexity of the interaction term, they are not regular like in preheating case. 



An exponential growth of the condensate for ever would be evidently catastrophic for this model. We 
will see later that when the radiation domination ends, the faster expansion of the Universe during 
matter domination stops the growth. On the other hand, if X has a short lifetime and the decay ends 
before the end of radiation domination epoch, production term becomes negligibly small and stop the 
growth. This is reflected on the approximate effective temperature of A particles obtained in Appendix 
1X1 The maximum amplitude of the condensate depends on the decay model, the lifetime of X, its 
coupling to <1>, masses, cosmological parameters such as Hq, and reheating temperature. We leave the 
numerical estimation of the quantities to a future work in which we solve the evolution equation of 
the condensate numerically. Nonetheless, classical treatment of the similar models [lOj shows that a 
condensate behaving like dark energy can be obtained for a large part of the parameter space without 
fine tuning. 

The self-interaction term in ()44p is highly non-linear: 

[G(^)>(x,y)]"— 1 - [G(^)<(x,y)]"-*-i (86) 

Note that we have separated the term corresponding to the classical potential. The term proportional 
to {(j)) is canceled out because the latter is null by definition. Assuming that the coupling A is small, 
we can linearize ()86p around xo = ^ which is the solution of (|74p when both couplings are zero. Even 
after linearization equation (j86p becomes a differentio-integral equation unless we replace all x's inside 
the integrand by xo^ except for i = term. Similar to the production term, the linearized interaction 
can be considered as a space-time dependent correction of the mass term which has the following 
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expression: 
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(87) 



Remind that according to (|3ip propagators of (j) have a term that depend on x and induces a dynamical 
mass for (p. This affects quantum corrections and their amphtude decreases with growing ip. We should 
also replace them with xo- The second complexity comes from the fact that it is non-linear in xo and 
its Fourier transform can not be described using - it depends on auto-correlation oiUkS. To simplify 
the expression we replace correlations with a simple multiplication. This is a rough approximation, 
but if Uk is dominant at /c ~ 0, this should not be a too bad approximation. Finally, after these 
simplifications the WKB correction term is modified to y\^'^"^ '■ 
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From preheating results [16] we know that even without taking into account quantum corrections, the 
first term in (|88p leads to a resonant amplification. However, the presence of the factor a^~" in front 
of the integral over rj in ()88p reduces the relative importance of this term at late times except when 
the self-coupling is much larger than the coupling of $ to X. But, in such models a strong clustering 
of the condensate is expected, which has not been observed in dark energy. Therefore, this type of 
models are not suitable as dark energy candidate. 

Quantum corrections of the self-interaction have only a minor effect on the total growth of the con- 
densate in this regime, because at lowest order they depend on rather than A in the classical 
term. We can make the same type of simplifying approximation that we have discussed above for 
the WKB integrals of production term and we obtain very similar expressions for the integrals - at 
least asymptotically. If we assume that the fraction of non-condensate ^> is negligible, T propagators 
(|119p and (jl20p contain only terms proportional to C and its conjugate. Equation (|118p shows that C 
also has the general form of the other expressions which are originated from ()55p and ()56p . Therefore 
after integration, the contribution of the quantum correction of self-interaction has the same general 
behaviour as the classical term and quantum corrections due to production - interaction with X. 
There is however, a major difference between production and self-interaction. The quantum field (/) 
has a dynamical mass that depends on the amplitude of the condensate. With the exponential growth 
of (f the effective mass of (j) increases and the amplitude of its propagator i.e. the cross-section of 
self-interaction decreases. This also affects {(pXA) in model (c) because its closed time path integrals 
contain a (p propagator. Therefore, in this regime the self-interaction can even decreases the amplitude 
of the condensate. It can be interpreted schematically as the recombination of (p and production of 
free cj) particles. 

We conclude that a classical scalar field - a condensate - can be formed during the radiation domination 
epoch with or without self-interaction. As the growth is exponential, the coupling constants must 
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be very small to prevent the over production of the condensate. Here we should remind that in 
the calculation presented here we didn't couple the expansion factor a(ry) to the evolution of the 
condensate and assumed that other constituent of the Universe, notably relativistic particles, dominant 
the energy density - the reason for calling this epoch radiation domination. In reality the evolution 
of the condensate and expansion factor are coupled and if the condensate dominates, the expansion 
factor grows more rapidly and reduces the growth of the condensate due to the presence of negative 
power of a{t) in the field equation and interactions. This fact becomes more clear in the next section 
where we study the evolution of the condensate during matter domination epoch. 

Before finishing this section we quickly review the behaviour of the two other models. Model (b) is 
very similar to model (a) but with an additional propagator in the closed time path integrals and 
an additional a^^ factor. Each one of them adds a factor of (??/%)~^^^ to integrals in w|''\ the 
analogue of W^°^ for this model. Therefore, W^''^ has 3 propagators in closed time path integrals and 
an additional {^/i]o)^^ factor. Although both these terms decrease with time, due to the exponential 
growth of Uk and Vk the condensate continues to grow exponentially but slower than model (a). 

In model (c) there are two expectation value due to the interaction of $ with X. The term {XA)c 
apparently looks like the similar term in model (a). But, the 99^ factor in ()24p makes it nonlinear. For 
obtaining an analytical solution, this term must be linearized and simplifications similar to what we 
explained for the self-interaction must be applied. The term {(j)XA) becomes linear if we replace the (p 
(or equivalently x) dependent terms in Gy with ipQ (xo)- As we mentioned above, the amplitude of this 
propagator decreases with the growth of (j). Therefore, the condensate evolution in this model should 
significantly deviate from the other models. In particular, we expect a stronger feedback between ip 
(equivalently x) and expansion factor a{rj). 



4.2 Matter domination epoch 



In this section we consider the matter dominated epoch. We limit ourselves to the time when the effect 
of dark energy is yet negligible. In the latter case we must consider the effect of (p{x) in the expansion 
of the Universe. This couples all the evolution equations and makes the problem insurmountably 
difficult. It is why we avoid this regime. 



The evolution equation of Xk for model (a) in this era takes the following form: 

d^kid^k2S^^\k-ki - h) 
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+ self-interaction = (90) 



This equation has an additional term due to a nonzero a^rj), and the expansion factor evolves according 
to (fml) in Appendix O It increases faster than radiation domination era, as f]'^ /iJq than rather than 
linearly. Similar to the previous section, we first neglect self-interaction. The solution of equation (|90|) 
without the interaction term is also explained in Appendix [Cj It has exact solutions for the special 
cases when m = or /c = 0. In Appendix O we use WKB approximation for the case of m 7^ and 
> 0. Because the interaction term in (|90p is proportional to x we can consider it as a time-dependent 
mass and use again the WKB technique to obtain an approximate solution: 
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As usual the integrals in ()92p are very complex and long. Therefore, in place of presenting all the 
details we only discuss the general form of their terms, and their asymptotic behaviour when r//?7o ^ 1- 

At late times we both A and X particles are non-relativistic and their temperatures can be considered 
to be zero. This makes the expression of propagators (|128p and (|129p much simpler. As we are only 
interested in the asymptotic behaviour of the condensate, we only consider the terms with the highest 
order of ?//?/o- Ii^ the integral over ?/ the dominant terms have the following forms: 



/ 
/ 



di—)\sma^ = ^-^^[e ""^"^(I-^ia^) - e'^e'^d, -ia\)] (93) 

,,1]' .7]'^ 7]'^ a~3 _Tjri 7 Vj'^ . lEi ,7 rj''^ , 
di—)—co^(^— = —^[(^ 7(0'*"— +e 6 7(-,-za— (94) 

with a = if3'j^ ± /3'^. Then the dominant terms in the integral over rj become: 

/T] 7 7]'^ rj'^ f V ^ v'^ v'^ 

d{—h{::,'ia—)sm{l3—), / d{—)-f{-,ia—)cos{l3—) (95) 
Vo 3 Vo Vo J Vo 3 ?7g ?7g 

where /3 = ±/3^ ± /3'x- These integrals can be calculated analytically when the 7 function is replaced 
by its asymptotic expansion ^{7/3,iar]'^ /ri^) r(7/3) — {iar}''^ /r]QY/'^exp{—iarj'^ /t^q). Finally the 
approximate expression for Xk at late times and without considering self-interaction is obtained as the 
following: 

• 9 /■ 
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4(2vr)3^7^l^ ' Vl' 

Y^CUk,xh{-\,ia^)] (97) 

where Ca and are proportional to the distributions of A and X particles in the same way as what 
was described in Sec. 14.11 Because in this regime both these particles are considered to be non- 
relativistic, equation ()8ip is applied as their distribution. At late times 7 functions in ()97p approach 
a constant and x dependent terms i.e terms containing f^^') and Z*-^^ decay very rapidly, as (?7o/^)^ 
for terms containing one /. Therefore, at late times Xkiv) is an oscillating function and its amplitude 
decreases as rjo/r] with time. Consequently, ipk decreases as r]Q/-rf and we conclude that in the matter 
domination epoch where the expansion of the Universe is faster, the production of $ in the decay of 
X alone is not enough to compensate the expansion, and the density of the condensate will decrease. 
Evidently, the validity of this conclusion depends on how precise are the approximation considered 
here. Note also that (j96p is written for small \k\ and is not valid for \k\ — )• 00. For the latter case the 
solution ()133p for m = must be used. It contains also only oscillatory terms. 

Equation (j90p becomes non-linear if we consider also the self-interaction. Therefore, to obtain an 
analytical solution we must linearize it by replacing the interaction term with its linearized version 
the equation ([57|) . and in 3^fc — ^ = 3^fc + -2^^) where Z'^ is the same expression as Zj. in ([5^ 
written for matter domination epoch. After taking integrals, Z'^^ will include terms very similar to what 
we obtained for 3^^. This means that the solution of complete linearized equation also contains only 
oscillating terms, and therefore decays as rjQ/rj^ or equivalently to/t. Nonetheless, for the reasons 
we explain now this approach is not realistic. In fact when self-interaction is added to equation 
([90|) . it becomes a non- linear differentio- integral equation, because the propagators of $ depend on 
the condensate. Therefore, only a solution without linearization can give a realistic answer to this 
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problem. This is possible only numerically and we leave it to a future work. In the rest of this section 
we estimation the late time behaviour of condensate qualitatively. 

To perform a non-linear analysis of equation (j90p with self-interaction we first neglect quantum correc- 
tions. This means that we only consider the first term of ()86p for which the minimum of the potential 
is at origin. If for simplicity we neglect also the production term the evolution equation becomes: 

^' + (A:2 + a'ml - \)x + Aa^-^x^'H^) = (98) 

Using a difference approximation for derivatives but without linearization, we find that although at 
the beginning x can grow, whatever the initial conditions, at late times it approaches to zero. This 
means that this equation lacks a tracking solution. Another way of checking the absence of a tracking 
solution is the application of the criteria F = y ' V jV'l > 1 suggested as the necessary condition for 
the existence of a tracking solution [36]. For equation (jOSp r = n(n — < 1 for n > 0. This is 

a well known result that only inverse power-law and inverse exponential potentials have a late time 
tracking solution [37j . 

The situation is different when we add quantum corrections. The coefficient C in (jll9p and ()120p 
which determines the amplitude of the quantum state of the condensate depends inversely on x- In 
another word, there is a back reaction from the formation of the condensate on the propagators of (j). 
To estimate the time evolution of the solutions of equation (j95p we use the approximate solution of 
field equation ()137p and determine the linearized self-interaction term ()87p . Similar to simplifications 
applied to the evolution equation in radiation domination case, we consider x as an external parameter 
and when we want to solve ()95p . we neglect its space dependence and identify fj with rj. The counting 
of X exponents in (j86p shows that due to the presence of negative power of x in (|118p some of the terms 
in the quantum correction, will be proportional to a negative power of x- This means that quantum 
corrections play the role of an inverse power-law potential with varying coefficients. As mentioned 
above, this type of potential is one of the well studied candidates for quintessence models and have a 
tracking solution [37] . The only difference here is the time dependence of coefficients. Therefore, we 
must find the conditions under which these coefficient are constant or vary slowly. 

It is easy to verify that the term with « = in (|86p includes the highest negative power of x- Therefore, 
it dominantly determines the tracking behaviour of the field. Fortunately, it has also the simplest 
integral in the sense that when we assume x as an independent variable and integrate over rj' , the 
integral contains a linear power of x- After taking the Fourier transform, Xk is integrated out and 
the integrand contains only the propagators. For i > terms the integrand includes also x factors 
and the evolution equations ()44p - ()46p become differentio-integral equation. For simplicity here we 
neglect these subdominant terms. After integration the evolution equation of the condensate in matter 
domination epoch becomes: 

a,/3 = i/3i, i = -(n-l),...,n-l (99) 

where dots indicates subdominant terms, including the production term. Therefore equation (j99p 
presents the dominant terms of the evolution equation for all three decay modes considered in this 
work. Ignoring the time-dependence of coefficients, it is evident that the effective potential in (j99p 
satisfies tracking condition explained above because it is a polynomial of inverse powers of x- But 
there is no known existence condition for tracking when coefficients are time dependent. The best 
guess is that if some of the coefficients approach a constant or vary slowly, the solution can be roughly 
a tracker. In fact we notice that due to negative powers of ry/jyo in ([99|) many of coefficients become 
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negligibly small asymptotically. By counting the order of ??/r/o terms, using the asymptotic expression 
of incomplete 7 function, we conclude terms satisfying the following conditions have slowest variation 
and are constant or growing when 77/% ^ 1: 

a = -2/3, 17-6n + 2i^0 (100) 

The first condition eliminates the oscillatory factors, and the second one is the order of rj/r]Q factors. 
As i ^ n — 1, this condition is satisfied only for n ^ 3. This is the only model in which the condensate 
does not decay quickly and a tracking solution does possibly exist. The case of n = 4 is also interesting 
because although apriori it does not have positive index term, it includes terms that decay only linearly 
with time, and therefore can lead to solutions in which the condensate decays with w > —1, but enough 
slowly to be consistent with data. We should also remind that renormalization of the complete theory 
in general induces anomalous dimensions for the fields that somehow modifies the exponents of the 
bare theory. In addition, although the production term decays with time, at intermediate epochs 
where matter is yet the dominant component its effect can be significant. 

Considering all these uncertainties, we conservatively conclude that for small n's such as n = 3 and 4 
the decay models considered here seem to produce a condensate that grows exponential in the radiation 
domination epoch and asymptotically evolves to a constant field in matter domination era. We note 
that these values for the self-interaction order are the only renormalizable polynomial potentials in 
4-dimension space-times. When the energy density of the condensate becomes dominant the study of 
its behaviour is more complex because the evolution equations of the condensate and expansion factor 
0(77) become strongly coupled. Nonetheless the classical treatment of the same type of models in |10| 
showed that the tracking behaviour persists in this regime, at least until the matter contribution in 
the expansion is not completely negligible. 

The amplitude of the effective interaction term in ()99p depends inversely on mass. Therefore fields 
with large mass produces smaller condensate. This is consistent with our initial assumptions that 
heavy scalars X and A don't condensate. Moreover, we note that ^ 1. This means that even 
for relatively large A the asymptotic value of (p can be small. This increase the initial probability of 
the formation of condensate and at same time reduces this amplitude at late times and reduces the 
fine-tuning of couplings. 

4.3 Initial conditions for the condensate 

To complete the study of the evolution of condensate we must fix the initial conditions. It is specially 
important for the asymptotic value of the condensate which determines its density today, and therefore 
constrains the mass and couplings of As we do not have an explicit expression for the solution of 
(|99p the initial conditions do not add any information to what we have obtained so far. Therefore, 
the discussion of this section is for the sake of completeness and would be useful future numerical 
simulation of this model. 

We assume that before the decay of X particles there is no condensate or more strongly there is no 
$ particle. In this case the initial condition for x is trivial, xiv = Vo) = 0. But x'iv = Vo) general 
should be non-zero and in fact positive. Its value presents the initial production rate of the condensate 
and is related to the initial density and the decay rate of X particles. 

It is expected that the Boltzmann equations which determine the evolution of distributions depend on 
if but not on its derivatives, see Appendix El Moreover, they are first order differential equations and 
each of them needs only one initial condition which can be chosen to be the initial distribution of the 
species. If we solve the evolution equation of ip along with the Boltzmann equations for all the species, 
this single set of inputs is enough to completely solve the Boltzmann and the condensate equation. In 
fact, once a solution of the Boltzmann equations as a functional of (p is found, by taking its derivative 
one can determine the initial value of ip' . This operation provides the complete initial conditions for 
the condensate's evolution equation. Here however, we do not solve Boltzmann equations in details. 
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Thus, we can not follow this procedure and must guess the initial value of 99' based on some physical 
properties. We remind that Boltzmann equations present the conservation of the flux of particles in the 
phase space. Thus we use the particle number conservation along with the additional assumption that 
initially, all the $ particles join the condensate immediately. We also neglect their spatial anisotropy. 
Under this conditions we obtain the following condition for the initial value of x'- 



Another way of determining integration constants is to fix the initial and final value of x- However, 
for the reasons explained above it will over-constrain the model when the full model is solved. For a 
tracking solution of ([UU]) the initial conditions don't affect the late behaviour of the condensate but 
they determine its asymptotic value. 

5 Discussion 

Although we have not yet observed any elementary scalar field, they are believed to play important 
roles in the foundation of fundamental forces and phenomena that have shaped our Universe. On the 
other hand, we have observed the composite scalar fields and their condensation in condense matter 
where they are responsible for various phenomena, such as symmetry breaking, mass acquisition of 
photons, quantization of flux tubes, formation of topological defects, and many other exotic behaviours 
of matter. From these findings we have learned that the condensate has usually a simple potential 
which can be related to the interactions in the original Lagrangian - at least qualitatively. For instance, 
in the case of Cooper pairs in superconductors, the presence of a tp^ potential can be schematically 
interpreted as an elastic scattering between electrons inside two Cooper pairs due to electromagnetic 
interaction. At low energies it is seen as a point-like scattering of 2 incoming scalar particles to two 
outgoing scalars - similar to self-interaction of a scalar field. 

The potential of the condensate of a fundamental scalar field should also trace back to the Lagrangian of 
the quantum field. Many particle physics applications of scalar condensates only consider the classical 
order which correspond to the potential in the Lagrangian. Because Lagrangians are usually local, the 
effect of the classical term is also local. However, as we have demonstrated in the previous sections, 
in some circumstances the effect of quantum corrections can be very crucial. In fact the properties 
of the condensate discussed in the previous sections, and in Appendices [X] and [B] are related to the 
non-local properties of quantum mechanics, and by extension quantum field theory. The descriptions 
of the quantum state of a condensate suggested in Ref. [2] and its generalization in Appendix |Bl 
include infinite number of entangled particles. Non-locality of quantum mechanics assures that these 
particles feel the presence of each others even at largest cosmological distances, and therefore behave 
collectively at large scales. At early times when the Universe is dense and the probability of scattering 
between particles is large, the local effect of the classical potential as well as the production of $ 
particles by the decay of X particles, which at lowest order is like a classical scattering, controls the 
amplitude of the condensate and the distribution of free (f) particles. But due to the expansion of the 
Universe, the cross-section of interaction and scattering becomes smaller at late times, and non-local 
effects become dominant. The very small coupling of <I> with other fields assures the stability of the 
condensate. It could be destroyed if interactions were strong enough to wash out the coherence of the 
condensate state, i.e. made particles to behave individually and semi-classically. 

Apriori it should be possible to design experiments or observations capable of distinguishing between 
a condensate of quantum origin and a fundamentally classical field as dark energy. If the latter case 
is true, dark energy must be due to a modification of the general relativity which is believed to be a 
classical theory at least for scales k < Mpianck- In this case the classical field would have a purely 
geometrical origin. On the other hand, if dark energy is produced by a condensate we expect to 
see some quantum effects. For instance to be able to observe its quantum excitation, similar to the 
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excitation of a Bose-Einstein condensate |41j . Evidently, due to the extremely small coupling of dark 
energy, the production of such excitations in the lab or their observation in cosmological environment is 
extremely difficult, if not impossible. Nonetheless, with the progress of our understandings in condense 
matter physics and related technological advance, there is hope that one day such an exploration 
become possible. 

If numerical simulations confirms our conclusions about the order of self-interaction potential, this 
would be a very significant result. In quantum field theories the dimension of an interaction term 
determines its renormalizability. For scalar fields in 4-dimension spacetimes with n ^ 4 are 
renormalizable. In fact except this physically motivated requirement there is no other rule to constrain 
the self-interaction in a quantum field theory. Interestingly enough, these values for n correspond 
exactly to models for which a late time cosmological condensate seems to excite. In the n = 3 models 
the coupling has a mass dimension of one which is considered to be the vacuum expectation value of 
another field, or to be proportional to the Planck mass, the only dimensional fundamental constant 
we know. In the latter case the field ^> is probably related to quantum gravity models. In addition, it 
must be a singlet or a 1-form in the group manifold of the symmetry group of the model, otherwise 
it breaks the symmetry. On other hand, in a n = 4 model the coupling constant is dimensionless 
and $ can be in a nontrivial self-conjugate representation of the symmetry group. These observations 
help to constrain the candidate particle physics models of dark energy. We can also put a rough 
lower limit on the self-coupling of <I>. r$ the interaction width of must be larger than Hq if <I> is 
not yet freezed-out. For n = 4 self-interaction this means A<j, > ^/Hqiti^, and for n = 3 potential 

A^/Mpianck > VHom^/Mpianck [E]- 

In summary, we used quantum field theory techniques to study the condensation of a scalar field 
produced during the decay of a much heavier particle in a cosmological context. Such a process had 
necessarily happened during the reheating of the Universe. It can also happen at later times if the 
remnants of the decay don't significantly perturb primordial nucleosynthesis. We showed that one 
of the necessary conditions for the formation of a condensate is its light mass. By considering three 
decay models and a power-law self-interaction potential of arbitrary positive order, we showed that 
the self-interaction has an important role in the cosmological evolution of the condensate and its 
contribution to dark energy. In particular, we showed that only a self-interaction of order 3 or 4 can 
produce a stable condensate in matter domination epoch. These results are obtained analytically and 
by considering a number of simplifying approximations. Therefore, they need confirmation by a more 
precise calculation, which in the face of complexity of this model, must be numerical. With little 
modification or adaptation, most of the formulation and results of this work are applicable to other 
cosmological phenomena which are based on a scalar condensate. 

We finish this section by reminding that if dark energy is the condensate of a scalar field, 
the importance of the quantum corrections in its formation and its behaviour found here 
is the proof of the reign of Quantum Mechanics and its rules at largest observable scales. 



Appendixes 

A Free field Green's function on non-vacuum states 

In canonical representation, a free scalar field $ can be decomposed to creation and annihilation 
operators on the Fock space: 

T = a(7y)$(x) = ^[Ukix)ak+UUx)al] , [ak,al] = 5kk' [ak^a^] = [4,4,] = (102) 

k 
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where Uk{x) = Uk{r])e ^^'^ is a solution of the free field equation (f5T]) . The quantization of $ imposes 
the following relation on Uk{x): 



U{x)W{y) -U{x)U*{y) = —6^^\x - y) 

a 



(103) 



A Fock state I*!') is constructed by multiple application of the creation operator a\ on the vacuum 
state |0): 



OfclO) = 0, VA; , \k1k2 ...kn) = a\^a\^ ■ • -^^lO) 

kik2...k„ 

The 2-point free Green's function of ^, it can be written as: 
iGF{x,y) = {^\T^x)^y)\^) = 

Yl ^kk^^k,k2..±S K{x)Uk{y)Q{xo-yQ)+Uk{x)Ul{y)e{yQ-xo) 

k i kik2...kn 



(104) 
(105) 



+ 



Uk{x)Ul{y)e{xQ - yo) +Ul{x)Uk{y)Q{yo - xq) 



(106) 



Prom ()106p we can extract the expression for advanced and retarded propagators: 

iG>{x,y) = {^^{xmym = YT. E 5kkS^k,k2...kSul{x)Uk{y) + 

k i kik2---kn 

E 1 + E E hkS^k,k2...kS Uk{x)Ul{y) 

k i k\k2...kn 

iG<{x,y) = {^mxMym = Y,Yl E ^kk^^k,k2...k,fmxp*kiy) + 



(107) 



i kik2...k„ 



E -^ + E E ^kki\'^klk2...kr, 

k i kik2...k„ 



(108) 



Therefore, for a free fields the Feynman propagator Gf{x, y) on a non-vacuum state can be written as 
a linear expansion with respect to the independent solutions of the free field equation V(k{x)V(^{y) and 
V(^{x)Uk{y)- The contribution of a non- vacuum state |^') appears in the coefficients of the expansion. 

When the entanglement and interaction between particles are negligible, the N-particle wave function 
1^*) can be written as a direct product of the 1-particle states. In this case, the projection coefficients 
are ~ f{k,x), where f{k,x) is the classical energy-momentum and space-time distribution of 
particles. In fact, if we project ^' to the coordinate space we can express as a functional of 
Wigner function |42j : 



1*1 



^'*(x)^'(y) = ^'*(x-h— )^'(x- — 



-9 



(2vr)^ 



x + y 



(fpP{x,p)e-'P-^ , x='^^ , X = x-y 

(109) 

In the classical limit the Wigner function P{x,p) approaches the classical distribution function f{p, x). 
As explained in Sec. \2.'6\ these distributions can be determined in a consistent way from the classical 
Boltzmann equations ()35p or their quantum extensions Kadanoff-Baym equations [HEH]. However, 
Boltzmann equations of interacting species are coupled to each others and to the evolution equation of 
the condensate. Thus, a complete solution can be obtained only by numerical calculation. Nonetheless, 
we need the distribution of particles to be able to to even an approximate solution for the condensate. 
For this reason we simplify the problem by assuming that f^'^\p,x) for i = X, A, (j) have thermal 
distribution and the effect of interactions can be included in the variation of their temperature. 
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The number density of X particles decreases by a factor of exp(— (t — toj/r) where to is an initial 
time. We assume that X particles are non-relativistic since their production i.e. Tx <C mx- Note 
that the decay is a non-thermal process. But when it happens slowly, the deviation from a thermal 
distribution is small and can be approximated by an effective time-dependent temperature. As X is 
non-relativistic its number density is approximately proportional to t'^'^ . This means that due to the 
decay Tx decreases by a factor proportional to power 2/3 of the decay term. We must also take into 
account the effect of the expansion of the Universe. Therefore, under these approximations: 

where D{x) is the growth factor of fluctuations, ks is the Boltzmann constant, and ksT^'^^ao) = 

2/3 

The distribution of A particles which are produced during the decay of X is also non-thermal. They 
are expected to be relativistic at early times if rriA ^ mx and non-relativistic at late times. Assuming 
that soon after they production they behave as a perfect fluid, the conservation equation gives their 
density, and thereby we find an effective temperature for them: 



^bTa 



2 - - _ , - 1 



3^-MaPx(x)^ - exp(-^)] - ^ + Ij 



4 



(111) 



where we have assumed a radiation dominated universe. When A particles become non-relativistic 
their effective temperature can be estimated as: 



^bTa 



2 {2tt)\ t-to 
-TjT-^ApxixJil - exp( )) 



(112) 



Note that in these approximations are obtained with the assumption that t ^ tjj, i.e. a very slowly 
decaying X. When t <^ tjj they should valid approximations when t <^ t. We emphasis again that 
the thermal distributions and temperatures calculated here are simple prescriptions when we can not 
solve all the evolution equations consistently. 

The case of ^ particles is somehow different because some of them join the condensate. Therefore, the 
effective temperature of non-condensate component depends on ip. Nonetheless, if we neglect their 
small coupling to other species, the total energy in the two components must be equal to the energy 
transformed to ^ during the decay of X. Their effective temperature have the same form as (jllip but 
with a branching factor that presents the fraction of energy transformed to non-condensated particles: 

A^^^A^^(l-^^^^±^) (113) 
Px 

where we have assumed that during X particles decay half of the energy is transferred to $ and the 
other half to A. We have also neglected the spacial fluctuations of the energy density of X and 

When we calculate the propagators of <j) we should take into account the contribution of all ^ particles 
in the wave function ^ , including the condensate. Therefore: 

|vl/(*)|2«/(*)(p,x) + /(^)(x) (114) 

where is the contribution of the condensate. Note that the separation of two components in 
(|114p is an approximation and ignores the quantum interference between them. It is valid if the 
self-interaction of ^ is weak and the non-condensate component decohere quickly. 

We ignore a general description for the wave function of the condensate component. Nevertheless 
special cases can be found, see e.g. [2] and Appendix |Bl These states are special cases in which 
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the addition of more particles does not change the state. In another word, they have a zero chemical 
potential. However, a condensate can be formed only in an interacting system, and in general the 
chemical potential is not zero. Using the description of the condensate wave function suggested in [2j 

1^^) ^ e-l^l'e^-^S |0) = e-l^l' ^(4)^0) 

It is easy to verify that this state satisfies the relation [2]: 

aol^'c) = Cl'^c) 

From decomposition of (p to creation and annihilation operators (jl02p we find: 

X{x) = a{^c\mc) = CUo{x) + C*U^{x) 



(115) 



(116) 



(117) 



Here we have adapted the original formula of [2] for a homogeneous FLRW cosmology. As x is a real 
field the argument of C is arbitrary and therefore we assume that C is real: 



C 



Uo{x)+U*^{x) 
X{x) 



(118) 



Assuming that the wave function of the condensate can be factorized, it is clear that in the expressions 
P07|) and (|108|) for the propagators, the condensate contributes only in terms in which at least some 
of the momentums are zero. Neglecting the interaction between particles, the wave function of these 
particles can be factorized and expressed as f^^\p,x) (see e.g. [3S])- The advanced and retarded 
propagators of (j) can be written aqj: 



1 



(2vr)= 



f^^\p,x) U;{x)Up{y) + (1 + /^<^^(p,x)) Up{x)U;{y) 



\C{x)\^U^{x)l(o{y) + {\C{x)\-' + l)Uo{x)U^{y) 



(2vr)= 



d^p 



/W(p,x) Up{x)U;{y) + (1 + /W(p,x)) K{x)Up{y) 



+ 



+ 



\C{x)\^Uo{x)U^{y) + (|C(x)|2 + l)U^ix)Uo{y) 



(119) 



(120) 



where E"^ = P"^ + m? . For fields X and A which do not have a condensate 11191 and 11201 can be use 
with C = 0. 

In Sec. m we used the Fourier Transform (FT) of the propagators in the quantum correction of the 
condensate evolution. To obtain the transformation of propagators in (|119p and ()120p we write them 
as the following (for C = 0): 



iG>{x,y) 



(fp 



fip,x) WJrjppir^yP-'' + (1 + fip,x)) Up{r])U;iv')e 



-ip.X 



(121) 



It is clear that G^{x,y) depends on both X and x, in contrast to Minkovski space-time where it 
dependens only on X. Nonetheless, it can be factorized to terms that depend only on one or the other 
coordinate. Therefore, the FT with respect to X is defined as: 



G>{x,X) 
G>{x,ri,r]') 



1 



(2^)2 



(fpG>{x,ri,T]')e' 



-ip.X 



f{-p,x) WJvP^pirj') + (1 + /(p,x)) U^mM) 



(122) 
(123) 



""it is evident that expression ()115|) is inspired from Bose-Einstein condensates in condense matter in which bosons 
with similar quantum numbers share the same energy state. Although this state satisfies the definition of a condensate 
according to no constraint on energy states or else can be concluded from the latter definition. Therefore definition 
(|115|l is a special case. 

®Note that for the condensate the momentum vector k presents its Fourier transform. By contrast, in /(fc, x) which 
is a classical distribution, k presents the momentum. 
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We remind that propagators are used to determine the closed time path integrals and thereby the 
expectation values in the evolution equation of the condensate. This equation is a partial differential 
equation, and we need to take its FT with respect to one of the variables x to solve it. If G{x, y) 
had a translation symmetry, the FT with respect to X variable was equal to the FT with respect 
to X or y up to a sign. But in an expanding non-empty universe there is no translation symmetry. 
Therefore, strictly speaking one can not use (jl23p in place of FT with respect to x. But, it is straight 
forward to show that the FT with respect to x mixes the modes in the quantum correction terms and 
makes any analytical solution of the condensate evolution equation very difficult. On the other hand, 
the fact that Gp {x,ri,7]') factorizes to components that depend only on one of the variables x and X 
suggests that as an approximation we can identify ()123p with the FT with respect to x, and treat x as 
an independent variable. Because the latter appears only in the energy distribution of particles, and 
U{x)U*{y) and its conjugate have translation symmetry this should be a good approximation. This 
allows to treat Gp {x,r],r]') like propagators in a Minkovski space-time. In Sec. [Hwe use Gp {x,r],r]') 
as the FT of propagators during determination of the closed time path integrals. 



B Generalized multi-condensate state 

Equation (jll5p which satisfies the definition ^ for a condensate (classical field) can be generalized 
in the following manner: Consider a system with a large number of scalar particles of the same type. 
The only discriminating observable is their momentum. The distribution of momentum is discrete if 
the system is put in a finite volume. Such setup can contain sub-systems similar to (jllSp consisting 
of particles with momentum k: 

l^fc) ^^fee^'=4|o) = Afc J]-^(4)^|0) (124) 

1=0 

where Ai. is a normalization constant. It is easy to verify that this state satisfies the relation: 

afc|^fc)7V = Ck\'i>k){N~i) (125) 

Therefore if — t- oo, the identity (|125p becomes similar to (|116p and the expectation value of the 
scalar field on this state is non-zero. However, if 7^ 0, the total energy of the system becomes 
infinite unless fc — )• and AA; — )• 0. Both here and the special condensate explained in Appendix [XI the 
mass of condensate is assumed to very small. In most realistic physical systems with a condensate, 
these conditions do exist. Therefore, we define a multi-condensate or generalized condensate state as 
a state in which every particle belongs to a sub-state of the form ()124p : 

N—>-oo 

\^GC) ^ Afce^'=4|0) = Y.A,Y^ -f (4r|0) (126) 

k k i=0 ' 

xix,^)^a{^Gc\mGc) = Y,CkUk{x) + C*M{x) (127) 

k 

The state \^gc) satisfies the equality ([125]). The coefficients Ck determine the relative amplitudes of 
the single-particle states with different momentum. Using (jl27p . the evolution equation of x determines 
how Cfc's evolve. 
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The extension of propagators ()119p and ()120p to this state is triviah 



E 



f^^\k,x)Ul{x)Uu{y) + (1 + f'^Khx)) Uk{x)Ul{y) + 



Ck{x)\' Ul{x)Uk{y) + (1 + \Ck{x)\') Uk{x)Ul{y) 



(128) 



iG'^<{x,y) 



E 



Uk{x)Ul{y) + (1 + Ul{x)Uk{y) + 



Uu{x)Ul{y) + (1 + Ul{x)Uu{y) 



(129) 



A simple example of such system can be the condensation of scalar particles in a potential well. Because 
only discrete energy levels are allowed, at equilibrium there can be a superposition of condensates with 
an effective mass (momentum) difference of Ame// = n/L where L is the size of the well. 

In the context of the model explained here the existence of a condensate in which particles have 
different energies is important because this means that $ particles do not need to lose completely their 
momentum to join the condensate. Such a state can potentially have applications in condense matter 
too, because in some sense it has simultaneously the properties of bosonic systems - condensation - and 
fermionic systems - a spectrum of energy levels. An example in the bulk excitation of a Bose-Einstein 
Condensate (BEC)-Bardeen-Cooper-Schreffer super fluidity [H]. 



C Propagators in matter dominated epoch 

In the matter dominated epoch the relation between comoving and conformal time is: 
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a 




ao 
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By applying ()13ip to the Green's function equation ()32p the field equation for the modes gets the 
following form: 



2 2 4 

m a^T] 

4 

Vo 







where Uk presents the solution for one of the field X or A. For two special cases of m = and k'^ 
this equation has exact analytical solutions: 



(132) 




vo ±2 V0-' 



/3' 



aor?om 
3 



For m 



■ 




(133) 



During matter domination epoch the masses of A and X particles are considered to be much larger than 
their kinetic energy and the Hubble constant. Moreover, in a cosmological context only large scales 
with k <^ m are under scrutiny. Therefore, we use A; = solution as the zero-order approximation 
and use the WKB-like techniques to find an approximation for k ^ case. Only the argument of the 
Bessel function in ()133p is mass-dependent. Thus, we replace it with a WKB-like integral: 



3 7?3 



aoVo 



J Vo Vo 



Vo 
. ^0 



1,2 

a? 

3^ " 
2m?r] 



(134) 
(135) 
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The Bessel functions J, i have analytical expressions: 




(136) 



After including ()136p to ()135p we obtain the following approximate expression for the solution of ()132p : 
Boundary and initial conditions explained in Sec. 13.11 are also applied to p37p . 

In the case of the field (j), the propagator (jl7|) includes also an additional mass term that depends 
on the condensate field (p, and therefore on rj. The effect of this term can be included in the same 
way as k'^ by replacing k'^'q^/m^rj terms with k'^tjo/m'^r] + A(n — l)y3'^~^/m|, in (jl37p . If at late times 
(fi — >• const - expected for a dark energy field - the contribution of this term does not vanish, in contrast 
to /c-dependent terms. Thus, it can be added to the mass of the field and the definition of /3'. This 
dynamical mass has the very important role of feedback on the growth of the condensate, see Sec. 14.21 
for details. 



D Propagators and condensate evolution in a fluctuating background 



Using the metric ()4ip . the Green's function equation ()3ip for the propagator Gpix^y) can be written 
as (we drop the field index for simplicity): 



(1-2^)2G'; 



il)'{G'j, - Gpa'/a) 



5'^ 



(1 - 2^) 



— - + (1 + 2^^) 2 didjGp 

(1 + 2^/^)2 
.4— n, ,n-2 



+ 



(1 - Aif')'2{a'^ml + {n- l)Aa^-"x""^) - (1 - 2^)^ — 

a 



Gf 



a 



Because (|138p is a linear differential equation and the metric fiuctuation ip is small Gp can be decom- 
posed to: 



GFix,y) = GUx,y) + AGFix,y) , 



AGF{x,y) 



G''Jx,y) 



< 1 



(139) 



where Gp{x,y) is the propagator in a homogeneous background i.e. the solution of equation (jl7|) . 
After inserting this definition in (|138p . we find the following equation for AGf- 

AG'; - 6'^d^djAGF + (a^ml + (n - l)Aa^-"x""' - — )AG'; = -i^S^x - y) + ({x - y) (140) 



C{x -y) = i;'{G'p - G%a'/a) + 5'^ {d^i^djG% + 2^d^djG%) - ip 



ih I 



a^ml + (n - l)Aa*-"x 



4— n, ,n— 2 



G% 
(141) 



The left hand side of (|14Up is similar to (|47p . thus both equations have the same homogeneous solution. 
If we neglect C, term, the only difference between (|140p and ()47p is the ip{y) factor in front of the delta- 
function in the right hand side of (|140p . Therefore, the solutions of these equations are the same up 
to a ip{y) factor. Because without <^ term the equation ()140p has the form of a Green's function, its 
solution in presence of a non- homogeneous term like Q is: 



AGf{x, y) = i^{y)G),{x, y) + ® C(^, v) 
The second term in (jl42p is of second order, and thus at first order in fiuctuations: 

GF{x,y) = {l + ^P)G^F{x,y) 



(142) 
(143) 
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Under the approximations considered in Sec. [Hwhen we calculate the closed time path integrals, the 
factors (1 + tp) of propagators can be included into A{k,x) in (|83p for radiation domination, and into 
C{k,x) in ()97p for matter domination. As explained in Sec. HI in this approximation the coordinate 
variable x is identified with x. 

The evolution equation for the classical component of $ must be also written for the background 
metric (|4ip : 



X 



i^'ix' - Xa'/a) 
(1 - 2^) 



T + {-, —J^^i^jX 



(1_ 4^2)5 1-2V' 
,4-n 



+ 



a^miil - i'4>^)2 -{l-2'4>) 



1 a 



+ 



n 



i=0 



n 



uj{x) 



'ag(XA) For ©-a 

g(XA2) = For ©-b 

^2(gx(XA)+g(TXA)) ForP-c 



(144) 
(145) 



In contrast to (jlSSp . for n > 2 this equation is non-linear and a decomposition similar to ()139p is not 
allowed. In the present work the aim of solving equation analytically obliged us in many places to 
neglect the coordinate dependence and non-linear terms. Therefore, in this work we neglect metric 
fluctuations in ()144p . This leads to equations 
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